Grain size dependence of barchan dune dynamics 
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Abstract. The dependence of the barchan dune dynamics on the size of the grains involved is investigated experimentally. 
Downsized barchan dune slices are observed in a narrow water flow tube. The relaxation time from an initial symmetric 
triangular heap towards an asymmetric shape attractor increases with dune mass and decreases with grain size. The dune 
velocity increases with grain size. In contrast, the velocity scaling and the shape of the barchan dune is independent of the 
size of the grains. 
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Our earth finds itself faced with the unpreventable cH- 
mate change [1]. The global warming makes the ice of 
the polar region melt and the deserts expand. To prevent 
further desertification, it is necessary to understand the 
dynamics of sand migration. The interaction between the 
moving sand and the driving air flow is highly compli- 
cated. Therefore it is desirable to investigate interactions 
of simple model fluid-grain systems. 

One of the fundamental dune systems is a barchan. 
Barchan dunes are crescent-shaped. They can be found as 
isolated dunes in unidirectional wind fields on bedrock, 
where sand is sparse [2]. These boundary conditions 
make their dynamics comprehensible for modeling, par- 
ticularly if the three-dimensional dune is cut in two- 
dimensional slices along the wind direction. Without 
losing the basic dynamics this simplification is done 
with so-called 'minimal models', which deal with two- 
dimensional barchan dunes [3, 4, 5, 6]. According to 
these models, the velocity scales reciprocal with the 
square root of the mass of the barchan dune [3] and re- 
ciprocal with the length of the windward side of the dune 
[4, 5]. Moreover, they predict the existence of a barchan 
dune attractor [7, 8], which has recently been demon- 
strated experimentally [9] . 

To test these minimal models it is necessary to 
study barchan dunes in well-controlled experiments in 
the laboratory. For a direct comparison with the two- 
dimensional models, it is useful to perform correspond- 
ing two-dimensional experiments. In former experiments 
we were able to verify the prediction of the velocity scal- 
ing laws of the minimal-models [10]. In those experi- 
ments the grain size of the dune's granulate has been 
kept fixed. In this paper, we address the question whether 
these results are still generally valid if the grain size is 
changed. 

A photo of our experimental setup is shown in Fig. 1. 




FIGURE 1. Photograph of the experimental setup. 



The main part consists of a closed flow tube machined 
from perspex which is filled with distilled water The 
height of the channel amounts to 60 mm, and its width is 
50 mm. The length of the straight section is 600 mm and 
the curves have an outer diameter of 500 mm. The water 
flow is driven by a propeller with a diameter of 45 mm, 
which is installed in the curve following the section of 
measurements. The flow direction is counter-clockwise 
on top view. 

We limit the section of measurements to approximate 
two-dimensional conditions, comparable to a Hele-Shaw 
cell [11]. For that purpose we use a black plastic insert, 
which constricts the channel to a width w = 6 mm. The 
aspect ratio of width to height amounts to 10. 

For the calculation of the Reynolds number Re of the 
6 mm wide gap we measure the flow velocity Vflo^ with 
an ultrasonic Doppler velocimeter {Signal Processing 
SA) and the kinematic viscosity v = 1 mm^/s of water 
at a temperature of 21 .7 ± 0.6C. The Reynolds number is 
kept constant for all measurements at Re = v^^^^w/v = 





FIGURE 2. Panel a) shows a side view snapshot of a fully 
developed barchan dune with mass m = 13 g. The glass beads 
appear white in front of a dark background. In panel b) the 
border of the dune is drawn as a black line. Five characteristic 
points of a barchan are shown as dots, and the length L and the 
height H are plotted as double arrows. 



27000. 

After the flume is filled with distilled water, glass 
beads are poured into the channel with a funnel at a fixed 
position upstream. The glass beads have a density of p = 
2.5g/cm^ and spherical shape. For the following exper- 
iments different grain size mixtures are used with diam- 
eters d ranging between 400-425 /xm, 450-500 ^m, 560- 
600 Aim, 710-750 Aim, and 800-850Aim. The selected 
range allows the comparison between the different grain 
sizes at a fixed Reynolds number At larger grain sizes 
Vflow has to be increased to trigger grain motion, and 
vice versa. The diameters are separated by multiple siev- 
ing. For the corresponding mean diameters d we obtain 
412 Aim, 475 /im, 580 Aim, 725 Aim, and 825 Aim, respec- 
tively. The different sizes of the barchan dunes are se- 
lected by weighing different amounts of glass beads. The 
selected masses m amount to 2. 17 g, 3.25 g, 6.5 g, 9.75 g, 
13 g, 16.25 g, and 19.5 g, with an error of ±0.005 g. 

A charge-coupled device (CCD) camera {Lumenera 
Lw 11059) with a horizontal resolution of 4008 x 2672 
pixels and a maximum frame rate of 5 fps is used for 
the measurements. The camera is placed in front of the 
straight part of the channel and records side views of 
the glass bead heaps and the developing barchan dunes, 
respectively. The lighting is done with a luminescent 
tube above the gap of the section of measurements to 
illuminate the upper surface. 

From pictures like Fig. 2a) the height profile of the 
barchan is extracted by finding the lowest gradient from 
dark to bright. In Fig. 2b) the border of the dune is drawn 
as a black line. We also indicate five dots which repre- 
sent the characteristic positions of a barchan: windward- 
bottom, crest, brink, lee-bottom and the center of mass 
(c.o.m.). The centroid of the two-dimensional border 
line indicates the center of mass of the barchan and 
serves to determine the mean dune velocity v. The other 
four points are determined by using the first and second 
derivative of the smoothed height profile, obtained by uti- 
Uzing a Gaussian filter The crest is located at the highest 




dune mass m (g) 

FIGURE 3. Final aspect ratio Hoo/Loo for different barchan 
sizes and different mean grain diameters d. The crosses denote 
the averaged final aspect ratio for each barchan mass m. The 
solid line is a guide to the eye. 



point of the dune. The brink lies at the maximum of the 
second derivative. The windward- and lee-bottom lie at 
the minima of the second derivative. 

The experiments always start with an elevated triangu- 
lar heap. During the time of measurements the heap be- 
comes lower and longer until its shape reaches its steady- 
state form. It has been shown that the migration veloc- 
ity of the dune takes a constant value instantaneously. 
Mass conservation is maintained all the time [10]. More- 
over, the relaxation time towards the attractor increases 
with mass and decreases with the flow velocity of the 
driving water flow. The final barchan dune attractor does 
not have the same shape for different flow velocities and 
different barchan dune sizes, that means that there is no 
scale invariance [3]. 

The final barchan dune attractor is reached when the 
ratio H/L becomes constant. The height H and the length 
L are defined as shown in Fig. 2b). To get a quantitative 
statement about the steady state, we fit the empirical 
function 

H //oo 

— = hoe ^ (1) 

L Loo 

to the temporal evolution of H/L [10]. The attractor is 
represented by the offset Hoo/Loo of Eq. (1). The aspect 
ratio a describes the deviation of the initial aspect ratio 
from its steady state value Hoo/Loo. The parameter T is the 
relaxation time for a triangular heap to achieve the stable 
shape of a steady-state barchan dune. 

The dependence of the final aspect ratio Hoo/Loo on 
the mean grain diameter J as a function of the size of 
the barchan dunes is shown in Fig. 3. The aspect ratio 
Hoo/Loo depends on the dune size. To get a better measure 
we averaged the final aspect ratios for each barchan size. 
The resulting mean aspect ratios are plotted as crosses 
in Fig. 3. For each grain mixture the deviation of Hoo/Loo 
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FIGURE 4. Relaxation time t for different barclian sizes 
and different mean grain diameters d. The crosses denote ttie 
averaged relaxation time fm for each barchan size. The solid 
line is a linear fit. 
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FIGURE 6. Relation between mean barchan velocity v and 
the mass of the barchan dime m for different mean grain diam- 
eters d. The lines are fits to v ~ m~". 
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FIGURE 5. Relaxation time T as a function of the mean grain 
diameter d for different barchan dune masses m. The crosses 
denote the averaged relaxation time Zj for each grain size. The 
solid line is a linear fit. 
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FIGURE 7. Mean barchan velocity v as a function of the 
mean grain diameter d for different barchan dune masses m. 
The crosses denote the averaged mean barchan velocity v* for 
each grain size. The solid line is a linear fit. 
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from the corresponding mean value seems to be random. 
So we suppose that there might be no coupling of the 
barchan shape, represented by Ho^/Loc, and the grain size. 

In Figure 4 the relaxation time T is plotted as a func- 
tion of the barchan size. To get a qualitative statement 
we average t for each dune size. For the seven differ- 
ent masses explored the resulting fm is shown as crosses 
in Fig. 4, and a linear fit is plotted to guide the eye. We 
suggest that for all different grain size mixtures the re- 
laxation time increases with the mass of the dune. 

To check the dependence of the relaxation time on the 
mean grain diameter we average t for each grain size 
mixture. The resulting fj and the original data is plotted 
as a function of d in Fig. 5. The mean relaxation time 
decreases with increasing mean grain diameter. 

In Fig. 6 the mean dune velocity is plotted versus 



the mass of the barchan dune. We fit our experimental 
data to the function v ^ m The resulting values for 
a are shown in Table 1. For the mean value we can 
assume a « 0.5 for all observed grain size mixtures. 
This corresponds to the theoretical prediction [3] and our 
former experimental result for a fixed grain size of d = 
560-600 nm and different flow velocities [10]. 

For the investigation of the dependence of the barchan 
velocity on the mean grain diameter we average v for 
each grain size mixture. The resulting v* and the origi- 
nal data is plotted as a function of d in Fig. 7. The mean 
barchan velocity increases with the mean grain diame- 
ter. This results from the fluidized surface layer whose 
thickness increases with the particle diameter c/ [12]. The 
larger the fluidized layer, the greater the amount of mass. 
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FIGURE 8. Relation between mean barchan velocity v" and 
the steady-state height of the barchan dune H„ for different 

mean grain diameters d. The lines are fits to v ~ HZ^ . 
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FIGURE 9. Relation between mean barchan velocity v and 
the steady-state length of the barchan dune Loo for different 
mean grain diameters d. The lines are fits to v ~ L^J. 



which is set in motion and shifts the center-of-mass. This 
also decreases the relaxation time as shown in Fig. 5. 

Figure 8 and Figure 9 show the relations between 
dune velocity and Hoo, respectively Loo, for different grain 



TABLE 1. Fit parameters obtained from the experi- 
mental data shown in Fig. 6, Fig. 8, and Fig. 9. 
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825 


0.52±0.01 


0.76 ±0.02 


1.00 ±0.02 


725 


0.54 ±0.02 


0.76 ±0.04 


1.04 ±0.04 


580 


0.55 ±0.03 


0.77 ±0.04 


1.02 ±0.07 


475 


0.47 ±0.03 


0.65 ±0.03 


1.08 ±0.05 


412 


0.53 ±0.03 


0.74 ±0.06 


1.10 ±0.04 



size mixtures. To get a quantitative relation we use the 
functions v ^ Hoo ^ and v ^ Loo^ as fits to the data points. 

Table 1 shows the results for the exponents j3 and 7. 
From these results we can assume that /3 = 0.74 and 
7= 1.05 within the errors of the measurement. This 
result corresponds to our former experiments with one 
fixed grain size and different flow velocities, on which 
the scaling laws v ^ //j" ** and ^ L^' are based [10]. 

To conclude, we investigated the dependence of the 
barchan dune shape attractor on its grain size mixture. 
We show experimentally that the relaxation time from 
a triangular heap towards the shape attractor decreases 
with the mean diameter of the grains and increases with 
the overall size of the dune. The dune velocity increases 
with the grain size. Moreover, the scaling laws for the 
velocity are independent of the grain size, as well as the 
shape of the attractor itself. 

The next step for understanding the barchan dune 
dynamics might be the investigation of the surrounding 
velocity field. Particularly the dependence of the surface 
roughness on the grain size and the coupling of surface 
and velocity field must be clarified. 
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